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Abstract. 

The extraction of the parameters of the supersymmetric Lagrangian from collider data is discussed. Particular emphasis 
is put on the rigorous treatment of experimental and theoretical errors. While the LHC can provide a first estimate of the 
parameters, the combination of LHC and ILC will be necessary to determine with high precision the parameters of the MSSM. 
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INTRODUCTION 

Supersymmetry with R-parity conservation is an attrac- 
tive extension of the standard model. It predicts a light 
Higgs boson mass in agreement with the electroweak 
precision data, and provides a candidate for dark mat- 
ter, the lightest neutralino. Moreover it provides the path 
to grand unification near the Planck scale. 

Among the various models presently being studied, 
two stand out: mSUGRA as an example of a model 
with few parameters defined at the GUT scale and the 
MSSM with many parameters which is defined at the 
electroweak scale. The supersymmetric particles are pro- 
duced in pairs, (cascade-) decaying to the lightest su- 
persymmetric particle (LSP), most naturally the lightest 
neutralino. The characteristic signature is missing trans- 
verse energy. 

In the following sections first the potential measure- 
ments at colliders, the LHC, the proton-proton collider 
with a center-of-mass energy of 14 TeV, and the ILC 
an e + e~ linear collider with a center-of-mass energy of 
up to 1 TeV, will be discussed. Then the reconstruction 
of the fundamental parameters will be described. The re- 
sults are summarized in the last section. 



MEASUREMENTS 

Supersymmetry can lead to a wealth of different signa- 
tures at colliders. To enable comparisons between ex- 
periments and/or colliders sets of parameters are defined 
to represent typical signatures. Most commonly known 
is the point SPSla [1] which has been studied in de- 
tail in [2]. Points with similar phenomenology have been 
studied in ATLAS (SU3) and CMS (LM1 Q). A sum- 
mary of the points is shown in Table [T] In the follow- 
ing SPSla will be studied as an example. The parame- 
ters of this point lead to squarks with masses of the or- 



TABLE 1. Summary of mSUGRA parame- 
ter sets with similar collider phenomenology. 
fl is positive for all cases. 





SPSla 


SPSla' 


SU3 


LM1 


mo 


100 


70 


100 


60 


my 2 
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der of 500 GeV /c 2 , light sleptons ( 1 50-200 Ge V /c 2 ) and 
the lightest Higgs boson mass being compatible with the 
LEP bound. 

An important aspect of collider phenomenology at the 
LHC are long decay chains such 
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The final state therefore contains at least a hard jet and 
two opposite sign same flavor leptons, where at the LHC 
lepton usually is restricted to be an electron or muon. In 
this decay chain five edges and thresholds can be cal- 
culated and reconstructed [2]. The expressions for these 
depend only on the four intervening masses. Therefore 
as the system is overdetermined the masses can be re- 
constructed, either via toy MC or fitting, without the use 
of assumptions on the underlying theory. Further signa- 
tures, e.g. the squark-R and the sbottoms, provide a total 
of 14 observable and measurable particles at the LHC. 
Typically the systematic error on measurements at the 
LHC coming from the jet energy scale is 1 % and 0.1% 
for the lepton energy scale. With integrated luminosities 
of up to 300 fb _1 the statistical error in many cases is 
smaller than the systematic error. 

The scenario of SPSla has been analyzed to great 
detail, providing to this day the benchmark point for 
studies at the ILC and LHC. More recent studies have 
increased the robustness of the results. The experiments 
have moved from fast simulations to full simulation and 



reconstruction This includes the detailed description of 
detectors as shown at this conference [0, 

While the determination of masses at the LHC can 
be performed most accurately by analyzing long decay 
chains, the ILC can measure, either via threshold scans 
or direct reconstruction the masses of essentially all kine- 
matically accesssible particles. This typically leads to a 
precision at the permil level. As the LSP mass is mea- 
sured more precisely at the ILC than at the LHC, one can 
insert the ILC measured LSP mass into the mass determi- 
nation at the LHC thus reducing the error on the squark 
mass measurement as shown in [71]. 

RECONSTRUCTING THE 
FUNDAMENTAL PARAMETERS 

Collider measurements are sensitive to several parame- 
ters and different observables depend on different combi- 
nations of the fundamental parameters. The system to be 
solved is therefore strongly correlated. A global Ansatz 
to make use of all information in an optimal way is there- 
fore necessary. 

Precise theoretical predictions are necessary in or- 
der to reconstruct the parameters. The ingredients [8] 
for such an endeavor, without attempting to be com- 
plete, are: precise mass calculations by SuSpect, SOFT- 
SUSY and SPheno (H [M [H ■ Branching ratios are pro- 
vided by SUSY-Hit which includes SDecay for the de- 
cay of supersymmetric particles and HDecay for the de- 
cays of the Higgs bosons lfl2l [l3l [Till . e + e~ cross sec- 
tions are calculated by PYTHIA and SPheno (H [Till , 
while NLO prot on-p roton cross sections are provided by 

ProspiM2.oEll3EII3l. 

To extract the parameters, the experimental data are 
analyzed by making use of these individual packages 
with a proper treatment of statistical, experimental sys- 
tematic and theoretical errors. Pioneering work on pa- 
rameter extraction with upward running to the GUT scale 
was done in 0. l2C}j. The Fittino 1I2TI1 group and the SFit- 
ter collaboration 12211 have worked on methods for the 
search for minima and the proper extraction of errors. 
More recently GFitter has presented a new package for 
the electroweak fit. The package Super-Bayes originated 
in the study of the dark matter aspect of supersymme- 
try JH. 

The standard model electroweak measurements with 
the addition of the WMAP [24] measurement of the relic 
density already allows to delimit interesting regions of 
parameter space without the direct observation of super- 
symmetric particles. In particular as shown in ll25l l26ll . 
the supersymmetric fit of mSUGRA results in a predic- 
tion for the lightest Higgs boson mass ofm h = 1 10+f d= 
3GeV/c 2 , pushing the Higgs boson mass closer, with 



TABLE 2. List of the best log-likelihood values over the 
mSUGRA parameter space using only the LHC measure- 
ments 1 22]. 
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respect to the standard model fit, to the limit of direct 
searches at LEP of 114.4 GeV/c 2 (zft. The distribution 
of the allowed cross sections for supersymmetric parti- 
cles at the LHC was shown in ll28tl . 



mSUGRA parameter determination 

Under the hypothesis that supersymmetric particles 
have been discovered and measured at the LHC and ILC 
and their discrete quantum numbers (e.g. spin) have been 
determined, there are two separate issues in the actual 
determination of parameters that need to be addressed: 

1. Can the correct SUSY parameter set be found? 

2. What is the precision of the determination of the 
parameters? 

mSUGRA provides a good testing ground for studying 
the techniques to answer these questions. The disadvan- 
tage of mSUGRA is that it starts with universal param- 
eters at the GUT scale, adds the RGE-extrapolation to 
the weak scale, and that it is not the most general La- 
grangian. 

To address the first question, SFitter performed 
300 toy experiments where the starting point of the 
simple MINUIT ll29ll fit was far away from the true 
values. The fits converged to the true values of SPSla 
already for the case where only LHC measurements, 
thus those with the largest errors, are available. 

Beyond the simple fit, powerful algorithms have been 
developed to sample multi-parameter space: The Fittino 
group has implemented simulated annealing which al- 
lows to escape out of secondary minima where a simple 
fit would have been confined to a wrong parameter sub- 
space. SFitter uses weighted Markov chains which allow 
for an efficient sampling of high dimensional parame- 
ter space. A full dimensional exclusive likelihood map is 
produced with the possibility of different kinds of pro- 
jection: marginalization (the Bayesian approach) or the 
profile likelihood (frequentist approach). Markov chains 
are also able to identify secondary minima. The comput- 
ing intensive method of dividing the parameters into a 
grid has also been implemented. 

The two-dimensional profile likelihood in the nuy- 
m\j2 plane is shown as an example in Figure Q] The 



TABLE 3. Best-fit results and errors for mSUGRA at the LHC (endpoints) and 
including ILC measurements taken from l22ll . 
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FIGURE 1. SFitter output for mSUGRA in SPSla. Two- 
dimensional profile likelihood x 2 over the triQ-m l n plane using 
only the LHC measurements. 



minimum at the correct nominal values is clearly visible. 
The search for secondary minima is illustrated in Table[2] 
Indeed, even in mSUGRA alternative minima can be 
found. It is interesting to note that the third minimum 
results from the interplay of the Aq parameter and the 
top quark mass (measurement error 1 GeV/c 2 ). This 
illustrates the necessity to take into account not only 
the supersymmetric measurements, but also the standard 
model parameters as part of the parameter sets. It is 
also a tale of caution: while the secondary minima can 
be discarded easily in this case by using the % 2 value, 
it is obvious that the unambiguous identification of the 
correct parameter set will become more complex in the 
MSSM case. 

To determine the correct errors on the parameters, 
not only experimental errors, but also the theoretical 
uncertainties, e.g. on the mass calculations, have to be 
taken into account. The RFit scheme (3011 : 
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(theo) 



K-rfi|<0) H (2) 
U-dilXT^ (3) 



defines the % 2 contribution as zero within the range of 
the theoretical error. This is appropriate as the theoretical 
errors are non-Gaussian. This treatment insures that no 



parameter value within the theory errors is privileged. 
Typical values for the theoretical precision are 3 GeV/c 2 
on the Higgs boson mass 13 ill . 1% on the masses of 
non-colored particles and 3% on the masses of strongly 
interacting particles lf22ll . 

Based on the mass determination at the LHC, the fun- 
damental parameters can be determined with a preci- 
sion at the percent level. Using the edges instead of the 
masses, thus using the experimental observables directly, 
the precision is improved by a factor 3 to 8. When cal- 
culating the masses from the edges, correlations among 
masses are introduced, as these are not available, the re- 
sults differ. A second remark is that using the correlation 
of the experimental systematic errors (energy scale) in- 
fluences the error at the level of 25% to 50%. Thus it will 
be important for the experiments to control these correla- 
tions with good precision. The best strategy to determine 
the parameters is to start from the experimental measure- 
ments and not intermediate quantities. 

Once the ILC becomes operational, the emphasis will 
turn even more to the precision determination of the pa- 
rameters. The result of the error determination for the 
LHC, ILC and their combination is shown in Table [3] 
In general the ILC improves the precision by almost an 
order of magnitude and the combination of the two ma- 
chines is more powerful than any single one. In partic- 
ular the measurement of tanj3 is improved as the neu- 
tralino/chargino sector is completely measured and the 
heavy Higgs bosons are within the kinematic reach of the 
ILC and can therefore be measured with high precision. 
Including theoretical errors does not change this picture. 
However it is interesting to note that they cannot even 
be neglected at the LHC alone. Therefore it is important 
already for the LHC to improve the precision of the theo- 
retical predictions. The SPA convention and project 03211 
provides a framework for this work. 

An interesting topic to illustrate that the determina- 
tion of parameters will be an iterative process (even for 
mSUGRA), is the neutralino enigma at the LHC. The 
LHC will be able to measure three neutralinos, therefore 
there is an ambiguity as to which ones: Xi> %® or 
Xv %2' %4 or %v ^3' Xa 1 - The last combination can be 
ruled out easily by the % 2 of the fit. However the first 
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FIGURE 2. Profile likelihood for the neutralino sector from 
S Fitter taken from [22]. 



two combinations both give a rather reasonable x ■ mo 
and mj/ 2 are correctly determined, only tan/3 and Ao are 
off, but they are not measured very precisely. However, 
the pre-determined parameters can be used to predict the 
mass and the branching ratio of the missing neutralino. 
As more are predicted than %®, the LHC alone can 
iteratively solve this problem. 



MSSM and extrapolation to the high scale 

The MSSM is more complex to analyze than 
mSUGRA as there are many more parameters to be 
determined. Typically 19 parameters are to be measured 
in the MSSM. In such a difficult environment not one 
single technique of finding the right parameter set but a 
mix is necessary. Taking the SFitter analysis as an exam- 
ple: a multi-step procedure of Markov chains alternating 
with MINUIT is used to analyze the MSSM. 

Three neutralino masses and no chargino masses are 
measured at the LHC in the point SPS 1 a. This results in a 
8-fold ambiguity in gaugino-Higgsino subspace. The % 2 
values in these points are essentially degenerate. Figure|2] 
shows the distribution of the inverse % 2 for M\. Four 
peaks are clearly seen. They correspond to the central 
values of the eight solutions. The solutions for positive 
and negative jA are too close to be distinguished in the 
Figure. 

When the ILC measurements are added, all parameters 
can be determined. Moreover, the precision of the mea- 
surement of the parameters is improved. The results of 
the error determination of the LHC and ILC individually 
as well the combination of the two is shown in Table [4] 
including flat theory errors. 

Having determined the MSSM parameters at the elec- 
troweak scale, the natural extension of the analysis (in 
contrast to mSUGRA) is the evolution of the parameters 
to the GUT scale as shown in lEojEUE^ . An example is 
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FIGURE 3. The upward renormalization group running of 
the inverse of the three gaugino masses in the MSSM as func- 
tion of the energy scale in GeV I7I I20I I2TI1. 



shown in Figure[3]for the inverse of the gaugino mass pa- 
rameters. As expected the three parameters unify similar 
to the unification of the gauge coupling constants at the 
GUT scale (10 16 GeV). This observation holds true also 
for the unification of the mass parameters in the sfermion 
sector as well as the tri-linear couplings. Thus by mea- 
suring the MSSM parameters at the electroweak scale, 
a window is opened to study the physics scenario at the 
Planck scale. 



SUSY with heavy scalars 

The point SPS la is a parameter set favorable for both 
the LHC and ILC. However even in difficult scenarios 
such as Split-SUSY HHEl (or SUSY with heavy 
scalars [36]) parameters of the underlying theory can be 
determined. 

The following scenario has been studied j3~7l [38ll : the 
scalar mass scale is set to at least 10 4 GeV/c 2 , thus the 
scalars are out of reach at colliders. The neutralinos and 
charginos have masses less than a TeV/c 2 . In the Higgs 
sector only the lightest Higgs boson is observable. 

The parameters to be determined are: ji the Higgs 
mass parameter, tan j3 the mixing angle in the Higgs sec- 
tor at the scalar mass scale, the gaugino mass parameters 
and the tri-linear coupling at the scalar mass scale. The 
scalar mass scale is fixed in this analysis. 

While the standard long decay chain is not available as 
the squarks are too heavy, the gaugino sector provides ad- 
ditional observables. The tri-lepton signal from the pro- 
duction and decay of the lightest chargino and the next- 
to-lightest neutralino provide information on the mass 
difference between the %® and %\ - Abundant production 
of gluinos is observed, and the cross section is used in 
the fit. 



TABLE 4. Results for the general MSSM parameter determination in 
SPSla assuming flat theory errors. As experimental measurements the 
kinematic endpoint measurements are used for the LHC column, and the 
mass measurements for the ILC column. In the LHC+ILC column these 
two measurements sets are combined. Shown are the nominal parameter 
values and the result after fits to the different data sets. All masses are 
given in GeV 1231 . 
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In the limit of infinite statistics, but taking into account 
theory errors, tan/3 is undetermined. By contrast, other 
parameters could be measured with a precision of several 
percent. The theoretical error on the gluino cross section 
was taken to be 30%. It is interesting to note that if the er- 
ror were negligible, the error on the gaugino mass would 
be improved by a factor 10. This shows the importance 
of precise theoretical calculations at the LHC. 



Connection to cosmology 

As demonstrated in [39] the connection between parti- 
cle physics and cosmology can be established at the LHC 
and ILC. First determining the parameters of the La- 
grangian via the collider data, the parameters can be used 
to calculate ( ||40l 5H SI) the relic density £2/i 2 . Shown 
in Figure [4] is the probability density of the relic den- 
sity for SPS 1 a' which differs from SPS 1 a only slightly as 
shown in Table Q] In this point mo and Aq were adjusted 
in order to lower the relic density which is too large to 
be compatible with the WMAP measurements ll24ll . The 
improvement from the combination LHC+ILC is clearly 
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FIGURE 4. Probability density of the relic density £lh 2 for 
SPS 1 a' for the LHC and the LHC+ILC (H] . 



visible. 

At SPSla in the mSUGRA scenario, neglecting the 



theory errors, a precision of 3% can be expected from 
the LHC alone, with the ILC improving the precision 
by an order of magnitude. The precision is comparable 
to that of the futur Planck experiment (2% H Hi ). 
The agreement between collider and cosmology mea- 
surements will be crucial to establish supersymmetry as 
the real solution to the dark matter enigma. 

CONCLUSIONS 

The LHC alone can provide a wealth of measurements 
for supersymmetry. The combination of LHC plus ILC is 
the optimal choice for the precision determination of all 
supersymmetric parameters at colliders, allowing a stable 
evolution of the theory near to the Planck scale. Once 
the parameters are determined it will be interesting to 
confront the prediction of the relic density from collider 
measurements with the measurements of WMAP and 
Planck. With the LHC starting this year, supersymmetry 
could be just around the corner. 

ACKNOWLEDGMENTS 

It was a pleasure to attend Supersymmetry 2008 in Seoul. 
I would like to thank Klaus Desch, Remi Lafaye, Tilman 
Plehn, Michael Rauch, Laurent Serin, Mathias Uhlen- 
brock, Peter Wienemann and Peter Zerwas for help in 
preparing talk and manuscript. The manuscript was writ- 
ten in the stimulating atmosphere of the workshop "LHC: 
beyond the standard model signals in a QCD environ- 
ment" at the Center of Physics in Aspen. 

REFERENCES 

1. B. C. Allanach, et al. (2002), hep-ph/0202233] 

2. G. Weiglein, et al., Phys. Rept. 426, 47-358 (2006), 
|hep-ph/041036"4l 

3. G. L. Bayatian, et al., J. Phys. G34, 995-1579 (2007). 

4. P. de Jong, these proceedings. (2008). 

5. P. Wienemann, these proceedings. (2008). 

6. C. Autermann, these proceedings. (2008). 

7. G. A. Blair, W. Porod, and P. M. Zerwas, Eu r. Phys. J. 
C27, 263-281 (2003), |hep^ph/021005"8] 

8. B. C. Allanach (2008). T0805 ■ 20881 

9. A. Djouadi, J.-L. Kneur, and G. Moultaka, Comput. Phys. 
Commun. 176, 426-455 (2007), |hep-ph/021133T) 

10. B. C. Allanach, Comput. Phys. Commun. 143, 305-331 
(2002), |hep-ph/0104145| 

11. W. Porod, Comput. Ph ys. Commun. 153, 275-315 (2003), 
|hep-ph/ 0301101] 

12. M. Muhlleitner, A. Djouadi, and Y. Mambrini, Comput. 
Phys. Commun. 168, 46-70 (2005), [hep-ph/ 3 1 1167] 

13. A. Djouadi, M. M. Muhlleitner, and M. Spira, Acta Phys. 
Polon. B38, 635-644 (2007), [hep-ph/ 60 9292] 



14. A. Djouadi, J. Kalinowski, and M. Spira, Comput. Phys. 
Commun. 108, 56-74 (1998), |hep-ph/ 97 04 448] 

15. T. Sjostrand, S. Mrenna, and P. Skands, JHEP 05, 026 
(2006), |hep-ph/ 6 3 1 7 5| 

16. W. Beenakker, R. Hopker, M. Spira, and P. M. Zerwas, 
Nucl. Phys. B492, 51-103 (1997), [hep"-ph/96104"90] 

17. W. Beenakker, M. Kramer, T. Plehn, M. Spira, 
and P. M. Zerwas, Nu cl. Phys. B515, 3-14 (1998), 
|hep-ph/971045T] 

18. W. Beenakker, et al., Phys. Rev. Lett. 83, 3780- 
3783 (19 99, Erratum-ibid. 100:029901,2008), 
|hep-ph/9906298| 

19. T. Plehn (1998), |hep-ph/ 9 8 93 1 9| 

20. B. C. Allanach, etal. (2004), |hep-ph/0403133| 

21. P. Bechtle, K. Desch, W. Porod, and P. Wienemann, Eur. 
Phys. J. C46, 533-544 (2006), |hep^ph/051100"6] 

22. R. Lafaye, T. Plehn, M. Rauch, and D. Zerwas, Eur. Phys. 
J. C54, 617-644 (2008), 0709. 39851 

23. R. R. de Austri, R. Trotta, and L. Roszkowski, JHEP 05, 
002 (2006), [hep-ph/ 602 028] 

24. D. N. Spergel, et al„ Astr ophys. J. Suppl. 170, 377 (2007), 
|astro-ph/06034"49] 

25. S. Heinemeyer, these proceedings. (2008). 

26. O. Buchmueller, et al., Phys. Lett. B657, 87-94 (2007), 
107 07 . 34471 

27. R. Barate, et al., Phy s. Lett. B565, 61-75 (2003), 
|hep-ex/0306033] 

28. B. C. Allanach, K. Cranmer, C. G. Lester, a nd A. M. 
Weber, JHEP 08, 023 (2007), 0705.04871 

29. F. James, and M. Roos, Comput. Phys. Commun. 10, 
343-367 (1975). 

30. A. Hocker, H. Lacker, S. Laplace, and F. Le Diberder, 
Eur. Phys. J. C21, 225-259 (2001), hep-ph/0104062| 

31. G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich, 
and G. Weiglein, Eur. Phys. J. C28, 133-143 (2003), 
|hep-ph/021202'0] 

32. J. A. Aguila r-Saavedra, et al, Eur. Phys. J. C46, 43-60 

(2006) , |hep-ph/0511344| 

33. N. Arkani-Hame d, and S. D imopoulos, JHEP 06, 073 
(2005), |hep-th/04 0515 9| 

34. G. F. Giudice, and A. Romanino, Nucl. Phys. 
B699, 65-89 (2004 Erratum-ibid.B706:65-89,2005), 
|hep-ph/040608"8] 

35. W. Kilian, T. Plehn, P. Richardson, and E. Schmidt, Eur. 
Phys. J. C39, 229-243 (2005), |hep^pri/040808"8] 

36. N. Bernal, A. Djouadi, and P. Slavich, JHEP 07, 016 

(2007) , [07TJ5.14 96I 

37. E. Turlay (2008), 08 5 .22721 

38. M. M. Nojiri, et al. (2008), [080 2 .3 67 21 

39. E. A. Baltz, M. Battaglia, M. E. Peskin, and T . Wizansky, 
Phys. Rev. D74, 103521 (2006), hep~-ph/ 602 187] 

40. P. Gondolo, et al., JCAP 0407, 008 (2004), 
|astro-ph/0406204| 

41. G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov 

(2008) , |0803 . 23601 

42. G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, 
Comput. Phys. Commun. 176, 367-382 (2007), 
|hep-ph/0607059"l 

43. J. Lesgourgues, S. Prunet, and D. Polarski, Mon. Not. Roy. 
Astron. Soc. 303, 45-49 (1999), [astro-ph/ 98 07 02 0| 

44. A. Balbi, C. Baccigalupi, F. Perrotta, S. Matarrese, 
and N. Vittorio, As trophys. J. 588, L5-L8 (2003), 
|astro-ph/0301192| 



